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Abstract 



There is now strong evidence that the main contribution to the 
cosmic energy density is not due to matter, but to another component 
with negative pressure. Its nature is still unknown: it could be the 
vacuum energy, manifesting itself as a positive cosmological constant 
with w = Pj p(? = — 1 (ACDM model), a spatially inhomogeneous and 
dynamically evolving form of energy with — 1 < ii; < (quintessence, 
QCDM model), or a dark energy component with w < such as 
a quantized free scalar field (VCDM model). After presenting sim- 
ple redshift-distance formulae, which are useful for comparing obser- 
vations with theoretical predictions without numerical integration, I 
discuss the behaviour of different tt; < models in the context of the 
Alcock-Paczyhski geometric test and the statistics of gravitational 
lensing. 

Keyw^ords cosmology: theory — cosmology: observations 
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1 Introduction 



The observations of distant Supernovae (Riess et al. 1998, Perlmutter et al. 1999, 
Riess et al. 2001) indicate that the universe is accelerating, as expected for a pos- 
itive cosmological constant; in this case it is possible to reconcile a relatively high 
Hubble constant with an old universe, and a low luminous+dark matter density 
with the flat geometry implied by the recent CMB observations (de Bernardis et 
al. 2000). The vacuum energy density can act as a cosmological constant, but 
its low value and its magnitude comparable to the present matter density repre- 
sent an unsolved problem (see e.g. Sahni & Starobinsky 2000, Carroll 2001, and 
references therein). An intriguing possibility is that the cosmic energy density 
is dominated by something acting as a variable cosmological "constant" (Ratra 
&; Peebles 1988); this additional component to the cosmic energy density, now 
known as "quintessence", would be slowly evolving with time and spatially inho- 
mogeneous, and could be due to a scalar field evolving in a potential coupled to 
matter through gravitation (Caldwell et al. 1998). Various authors (see Efstathiou 
1999, Wang et al. 2000, Balbi et al. 2001) have shown that a flat model including 
a dominant component of quintessence (with equation of state wq ^ —0.65, see 
below) is consistent with observations, while Amendola (2000) has put constraints 
on a quintessence-dark matter coupling. 

The observational constraints are still quite uncertain and it is not possible 
to discriminate between the different model parameters (Podariu & Ratra, 1999). 
Moreover, Caldwell (1999) suggested that even models with w < —1 satisfy all 
the observational constraints; he named this hypothetical component "phantom 
energy" , as in his model there is a negative kinetic term in the Lagrangian of the 
scalar field (see also Schulz & White 2001). Nevertheless a w < —1 equation of 
state can also be obtained assuming that the vacuum energy is due to a quantized 
free scalar field of low mass (VCDM model, Parker & Raval 2001 and references 
therein). 

Changing the equation of state of the dominant energy component of the uni- 
verse affects the distance-redshift relation and all the distance-dependent physical 
quantities. The effects become important in deep redshift surveys. Unfortunately, 
in most cases it is not possible to find exact analytical expressions for the distance- 
redshift relation (e.g. McVittie G.C. 1965). Models with a cosmological constant 
are now well studied, and distances and other useful quantities have been obtained 
through numerical integration (Refsdal, Stabell Sz de Lange 1967) or using elliptic 
integrals (Feige 1992). An excellent approximation for A > flat cosmologies has 
been found by Pen (1999), while Kayser, Helbig & Schramm (1997) describe a 
method to calculate cosmological distances also for inhomogeneous cosmological 
models. Recently, Hamilton (2001) has also found formulae for the linear growth 
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factor and its logarithmic derivative. 

On the other hand, the alternative models have not been studied extensively 
(but after the submission of this paper, the situation has continued to improve; see 
e.g. Giovi et al. 2000). As simple expressions for the distance-redshift relation are 
always useful to perform a fast analysis, in this paper I present simple formulae 
for quintessence cosmological models. 

It is assumed that a) the Universe is homogeneous on large-scales, so that 
Priedmann-Lemaitre equations can be used -in the case of gravitational lensing 
filled beam (i.e. standard) distances will be adopted-, and b) the equation of state 
is constant. 

In section 2 the basic equations are briefly described, and some exact analytical 
formulae for the redshift-distance relation are explicitely shown in the trivial cases 
of one component with w < 0; it is also shown a formula with elliptic integrals 
for a realistic matter + quintessence cosmological model with w = —2/3, which 
satisfies all the presently available observational constraints. 

In section 3 and 4 I discuss the differences among different models with a dark 
energy component which are expected when applying respectively the Alcock- 
Paczyhski test and the statistics of gravitational lensing. The main conclusions 
are in section 5. 



2 Redshift-distance relations for models with 

w <0 

2.1 Basic equations 



For sake of clarity, the basic equations for calculating distances in relativistic 
cosmology are briefly resumed (e.g. Peebles 1993; Coles &: Lucchin 1996, Pea- 
cock 1999; see also Hogg 1999). Assuming that the universe is homogeneous and 
isotropic universe, we have the Robertson-Walker metric: 



where Sk{r) = r for k = 0, Sk{r) = sin(r) for = 1 and Sk{r) = sinh(r) for 
A; = — 1, and from Einstein's field equations the Priedmann equations are obtained: 



ds^ = c^dt^ - R^{t) dr^ + Sk{r){de^ + sin^ 9d<p^) 



(1) 




(2) 



R 
R 



^2^Pi(l + 3u'i) 



(3) 
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where pi is the energy density of component i and Wi = Pi/ pi(? is the corre- 
sponding equation of state. 

Defining the ratios = Pi/ Pc, where p^ = 3H'^/8ttG and 17^ = —k(P'/{HKf', 
we can rewrite equation @) as: Vltot = J2i = 1 — ^k- 

Moreover, dividing equation (^) by equation (^) we obtain the general expres- 
sion for Bandage's deceleration parameter q = —RR/E? = ^(1 — ^fc) + | Z]j ^iWi- 

The comoving coordinate r can be written as a function of 

dz' 



where: 



E{z) = ^nkii + zy + Y. niii + z)3(i+-.) (5) 

In what follows I will consider only models with one w < component and 
matter {w = 0), while radiation will be neglected. In the following, whenever 
useful, a subscript will indicate the range of values assumed by the equation of 
state: wq for values in the range —l<w<0, and wy for values less than —1. 

The comoving distance]] is therefore given by: 

= RoSkir) = ^-1—Sk{r) (6) 

where, as usual, the subscript indicates the value at the present epoch. The 
luminosity distance di is simply given by = dc{l+z), while the angular distance 
is dAiz) = del {I + z). 

The lookback time is also given by an analogous integral: 

I 1 
^^'^^IfoJo il + z')Eiz') 

Numerical integration (see Press et al. 1992) has been used to check the results 
obtained with analytical formulae and more generally when no analytical solution 
was available. 



2.2 One— component models 

Simple analytical expressions can be found only in special cases. Calculations are 
clearly simplified assuming w constant (a reasonable approximation as far as the 
equation of state changes slowly with time). 

^dc{z) corresponds to the transverse comoving distance Dm as defined by Hogg (1999). 
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For a flat universe (0^ = 0) with only one component, analytical solutions to 
the integral in equation (^) are very simple. Defining a = (1 + 3i(;)/2, we obtain: 



4(z) 



+ for w + -\l?, 



for w 



-1/3 



(8) 



Other simple solutions can be found by fixing the value of w. The values w = 
— 1/3 and w = —2/3 are of special interest because they correspond respectively 
to a frustrated network of cosmic strings and of domain walls (Bucher and Spergel 
1999); morevoer, the value ti; = — 2/3 is consistent with present observations. 

For w = —1/3 we have: 



sinh(^l-!^Qln(l+z)) 



Ho 



for <Q,Q <1 
for f^Q = 1 



and for w 



-2/3: 



(9) 



[2,1(^1+^-1) nQ = l 

(10) 

Figure la shows the distance-redshift relation for flat cosmological models with 
one component having different equations of state, and assuming Hq = 65 km/s 
Mpc~^. It is clear that a more negative equation of state implies a larger distance 
at a given redshift, and that the difference increases with redshift. 

2.3 A formula for the "standard" quintessence model 

Numerical integration of equation (Q) is necessary in the general -and realistic- 
case of a two-component model including both matter and dark energy. In the 
case of flat cosmological models, the comoving distance of a galaxy at a redshift z 
can be written using the hyper geometric function 2-^1 (oj b; c; z): 



dc{z) 



1 



--{ 2F1 



2F1 



X 

1 1 



1 1 

6^' 2' 



1 

Qw ' 



1 



X 



l-Qx 

3w 



{l + Z 



} 



(11) 



where Qx refers to the dark energy contribution. The above relation is con- 
sistent with the formula found by Torres and Waga (1996) for angular diameter 
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Figure 1: a) Comoving distances for flat cosmo logical models with one com- 
ponent. Dotted dashed line: w = —3/2; sohd curve: w = —1; dashed curve: 
w = —2/3; dotted curve: w = —1/3; long dashed curve: Einstein-de Sitter 
model, b) Comoving distances for an Qm = 0.3, open model (dotted curve), 
the Einstein-de Sitter model (long-dashed curve), and 3 fiat models, with 
= 0.3 and an additional negative pressure component Qx = 0.7 with the 
following equation of state: dashed-dotted curve: w = —3/2; solid curve: 
w = —1; dashed curve: w = —2/3. 
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distances. Assuming w = —2/3 it is possible to further simplify the above ex- 
pression, using only the incomplete elliptic integral of the first kind. A simpler 
expression can be of practical utility as this equation of state is consistent with 
the restricted range permitted by present observations (the best parameters for a 
flat universe are 0,^ = 0.33 it 0.05 and wq = 0.65 ± 0.07 according to Wang et al. 
1999). 

The formula is the following: 



dc{z) = < 



where: 



C 
C 

C 



2^(1' 71 



«(5(^),^) 
n(<7(0),^) 

(<,(0),^)-n(<,(z),^ 



< < 0.5 

0.5 <VLq<1, z< 

0.5 <Q.Q <l, z> 



and 



u{g: 



C 



9iz) 



71' 



1 



+ 1 + Z 



I-Qq 



-.dt 



(12) 



(13) 
(14) 



(15) 



V(l-t^)(l-t^/2) 

is the incomplete elliptic integral of the first kind, which can be easily com- 
puted with one of the standard routines available in many mathematical packages]^. 
u(l,l/\/2) = -f^(l/\/2) corresponds to the complete elliptic integral of the first 
kind, and u((/(0), l/\/2) depends only on Oq: when Qq = 0.5, u{g{0),l/ V'^) = 
K{l/V2). 

Figure lb shows dc{z) for the Einstein-de Sitter model, an open model with 
r^M = 0.3, and other 3 flat models with Qm = 0.3 and a negative pressure com- 
ponent with w = -2/3 (QCDM), w = -1 (ACDM), and w = -3/2 (VCDM). 
A synthesis of their main properties is given in table (1). A comparison between 
figures la and lb shows clearly that the difference between models becomes small 
in the presence of matter. For example, up to z ~ 3 the curve of the open model 
is very similar to that of the quintessence model with w = —2/3. 
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12 3 4 

Z 



Figure 2: a) Angular diameter corresponding to 1 Mpc (for Ho = 65 
km/s/Mpc) as a function of redshift for the same models as in figure lb. 
Dotted curve: = 0.3, open model; long-dashed curve: Einstein-de Sitter 
model; 3 flat models, with Qm = 0.3 and an additional negative pressure 
component fix = 0.7 with the following equation of state: dashed-dotted 
curve: w = —3/2; solid curve: w = —1; dashed curve: w = —2/3. 
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Table 1: Properties of different cosmological models 





(1,0,0) 


(0.3,0,0) 


(0.3,0.7,-2/3) 


(0.3,0.7,-1) 


(0.3,0.7,-3/2) 




0.5 


0.15 


-0.2 


-0.55 


-1.07 


Age (lO^het yrs) 


10.3 


12.2 


13.6 


14.5 


15.4 




1.25 


1.90 


1.59 


1.61 


1.55 






1.02 


3.29 


1.48 


0.79 




1.00 


1.06 


1.24 


1.23 


1.22 



Another interesting quantity is the angular diameter. The angular diameter of 
an object with physical diameter D \s 9 = ^^'^^^^ = ^^^{z) ■ Figure 2 shows the 
angular diameters corresponding to 1 K^^ Mpc for the same 4 models as in figure 
lb. Also in this case the differences between models with negative pressure are 
not very large, and it is apparent the similarity between the open and standard 
quintessence models. A detailed discussion of the angular size in the context of 
quintessence models can be found in Lima &; Alcaniz (2000). 

3 The Alcock— Paczynski test 

Alcock k, Paczynski (1979) showed that, at least in principle, it is possible to de- 
tect the presence of the cosmological constant A in a redshift survey through a 
geometric test which is independent of galaxy evolution. This method is based 
on the fact that the spatial distribution of galaxies obtained from the distance- 
redshift relations valid for an Einstein-de Sitter model is significantly distorted if 
the universe has a A component. Therefore a measure of the apparent anisotropy 
in known isotropic structures would give us the value of the cosmological constant. 
Phillips (1994) applied this method to the orientation of quasar pairs, without con- 
clusive results. His work has been extended to the quasar correlation function by 
Popowski et al. (1998). The main problem is that redshift distortions are also in- 
duced by peculiar velocities -/^-distortion (/? = $7^/6) on linear scales and fingers 
of God on small scales-, but the method might be successfully applied to future 
redshift surveys, with the simultaneous extraction of A and (3 from anisotropic 
power-spectrum data, as discussed by BaUinger, Peacock k, Heavens (1996). 

^For example EllipticF in Maple V or ellfin Numerical Recipes: notice that e/y requires 
the quantity asm{g) instead of g. 
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Z Z 

Figure 3: Flattening factor F{z) for different flat models; in figures a-c, the 
value of flM varies from 0.1 (top curve) to Qm = 0.9 (bottom curve) in steps 
of 0.1. a) QCDM (fi^ + VIq = 1, wq = -2/3); b) ACDM (^m + = 1); 
c) VCDM {VIm + nv = I, wv = -3/2); d) 4 flat models with = 0.3. 
Dashed curve: QCDM model with wq = —2/3; solid curve: ACDM model; 
dashed-dotted curve: VCDM model with wy = —3/2; dotted curve: VCDM 
model with wy = —3. 
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The distortion can be quantified by a flattening factor F{z) relative to the EdS 
model. In our notation, it is defined as: 



F > 1 indicates a flattening along the line of sight (see Ballinger et al. 1996). 
Only in special cases F{z) can be expressed analytically: for example, for a flat 
quintessence-only model (Oq = 1) and wq = —2/3, we have F{z) = \/l + z, to 
be compared to the = 1 case, with F{z) = 2(1 + z)[y/ (1 + z) — l]/z. 

Note that the formula valid at first order in z given by Ballinger et al. (1996) 
can be generalized to quintessence models: 

F{z) = i + ^[i-nM-{i + 3wQ)nQ]z + 0(^2) (17) 

Therefore, when wq = —2/3, the test is a direct measure of the difference 

I — I • 

In figure ^-c it is shown F{z) for flat models with w = —2/3, w = —1, and 
w = —3/2, while figure ^d shows 4 flat models with ^Im = 0.3. 

There are three main systematic effects due to the different equation of state 
which appear clearly in figure ^d. 

With a more negative w we have: 

• a) the amplitude of the maximum anisotropy becomes (slightly) smaller; 

• b) the maximum of F{z) shifts towards smaller redshifts; 

• c) after the maximum, F{z) decays faster. 

With a realistic fraction of matter the amplitude of the geometric distortion is 
not very large, as it amounts to about 20%. After the first analysis of Boomerang 
results, McGaugh (2000) stressed that a universe with ~ 1 and 0,m = would 
be consistent with results from the observations of distant type la supernovae and 
with the small amplitude of the second peak in the power spectrum of anisotropics 
in the CMB, if instead of CDM one adopted MOND (Milgrom 1983). In this 
case, we would expect a larger geometrical distortion, about 30% at z = 0.5 and 
50% at z = 1. However, the more extended analysis of Boomerang results, which 
has recently identified the first three peaks in the CMB angular power spectrum, 
shows that a standard flat model with = 0.7 and = 0.3 is consistent with 
primordial nucleosynthesis (Netterfield et al. 2001). 

An interesting aspect of the Alcock-Paczyhski test applied to the VCDM mod- 
els is that the maximum distortion is shifted towards relatively low redshifts: 
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Zmax ~ 0.32 if wy = —3, and even Zmax ~ 0.08 in an extreme case wy = —10. 
This well-defined maximum at moderate redshifts could be detected in on-going 
galaxy redshift surveys, or at least it should be possible to fix a lower limit to the 
equation of state, depending on our ability to disentangle the geometric distortion 
from the redshift distortion due to peculiar velocities. 



4 Statistics of gravitational lensing 

Gravitational lensing represents an important test of cosmological models, as it is 
sensitive to the presence of a cosmological constant (Fukugita et al. 1992) or to 
an equivalent component with negative pressure. Here we simply aim to examine 
the relative trend for the different models, without going in much detail (among 
the more recent studies see e.g. Kochanek 1996, Chiba & Yoshii 1999). 

The differential probability dr of a line of sight intersecting a lensing galaxy 
(modelled as a Singular Isothermal Sphere) at redshift zl in the redshift interval 
dzL is: 

E{zl) \hJ dAi^^zs) ^ ' 

For flat models, we have the analytical relation: 



r(.5) = f^dA^sY [j-J (19) 

where / measures the effectiveness of the lens in producing double images 
(Turner, Ostriker & Gott 1984); here its value is estimated assuming a Schechter 
luminosity function: 

/ = ^**''"r(« + - + i) (20) 

Figure ^ shows the optical depth as a function of the redshift of the source, 
calculated using the angular diameter distance. With the above relations, it is 
possible to compare the predictions of gravitational lensing for various cosmological 
models to observations. Torres & Waga (1996) have examined the behaviour of a 
model with w = —1/3. Here their analysis is extended to include the more popular 
value w = —2/3 for quintessence, and the case w = —3/2 for VCDM. I have used 
the same catalog, the HST Snapshot Survey (see Maoz et al. 1993), which has 
detected a number between 3 and 6 of gravitational lensing candidates out of 502 
observed quasars. Magnification bias has been taken into account as in Maoz & 
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Figure 4: Optical depth for different flat models with ^Im = 0.3. From top to 
bottom: dashed-dotted curve: VCDM model with wy = —3/2; solid curve: 
A model {wa = —1); dashed curve: quintessence model with wq = —2/3; 
long-dashed curve: Einstein-de Sitter model. 
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Table 2: Predicted number of lensed quasars and corresponding probabilities. 



(1,0,0) 


(0.3,0.7,-1) 


(0.3,0.7,-2/3) 


(0.3,0.7,-3/2) 


(0.5,0.5,-3/2) 


Lens no. 3.7 


10.8 


8.6 


13.5 


8.0 


Prob. 0.193 


0.012 


0.041 


0.002 


0.058 



Rix (1993), and I have adopted ^* = 0.015 h^Mpc-^, « = -1.2, a* = 200 km/s 

and 7 = 4. 

Table 2 resumes the main results, assuming 4 genuine lensing events in the 
HSTSS. It is well known that A > models predict a too large number of lensed 
quasars, as the optical length increases with | u; | (as shown in figure 4; see also 
Zhu 2000). While the wq = —1/3 quintessence model considered by Torres &: 
Waga (1996) could still be consistent with the data, the same is not true for the 
Wq = —2/3 model, and VCDM models perform still worse than the A model. For 
VCDM models, more acceptable results can be obtained if the relative contribution 
of matter is increased. For example, in a fine-tuned VCDM universe with equal 
contributions of matter and dark energy {^m = 0.5 and Jly = 0.5) the number of 
lenses becomes comparable to the quintessence wq = —2/3 model (see table 2). A 
large matter density results also from the analysis of radio-selected gravitational 
lenses and optically selected lensed quasars, which gives CIm > 0.38 assuming a flat 
universe and a nonzero cosmological constant (Falco, Kochanck &: Mufioz 1998). 

New surveys should clarify the reasons of this cosmic "discordance" with re- 
spect to the increasing observational evidence that the energy density of the uni- 
verse is dominated by a component with negative pressure: it could be due to 
systematic errors and uncertainties in the modeling of gravitational lenses (Cheng 
& Krauss 2000). 

5 Conclusions 

In this paper various aspects of negative-pressure cosmological models have been 

analysed: 

a) simple redshift-distance formulae (the equivalent of Mattig formulae) have 
been derived for quintessence cosmological models; 

b) it has been shown that the Alcock-Paczyhski test is particularly sensitive to 
w < —1 models, where the distortion has a maximum at a relatively low redshift; 
in the case of a quintessence model with w = —2/3, this test measures directly the 
difference | flq — | • 
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c) as the more negative is w, the larger is the optical depth, a quintessence 
cosmological model with w = —2/3 predicts a lower number of lensed quasars -^and 
is therefore more consistent with observations- than the corresponding A model, 
while the w = —3/2 model gives an even larger number of expected lenses, which 
can be reduced by increasing Qm- 

A more refined analysis is obviously required to test models with negative 
pressure and their consistency with the full set of observational data (for exam- 
ple, adopting tracking models with a variable equation of state, see Benabed & 
Bernardeau 2001). 

Oncoming deep redshift surveys, such as the VLT-VIRMOS Deep Survey (Le 
Fevre et al. 2001) and DEEP (Davis & Newman 2001), will sample distances larger 
than z ~ 1, so that it should be possible to obtain constraints to the cosmological 
parameters which will be complementary to other measurements. 
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